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ABSTRACT 

We make a detailed measurement of the metal abundance profiles and metal abundance ratios of the inner core 
of M87/Virgo observed by XMM-Newton during the PV phase. We use multi temperature models for the inner 
regions and we compare the plasma codes APEC and MEKAL. We confirm the strong heavy elements gradient 
previously found by ASCA and BeppoSAX, but also find a significant increase in light elements, in particular O. 
This fact together with the constant O/Fe ratio in the inner 9 arcmin indicates an enhancement of contribution in 
the core of the cluster not only by SNIa but also by SNII. 

Subject headings: X-rays: galaxies — Galaxies: clusters — Galaxies: individual: M87 — Galaxies: abundances 



1. INTRODUCTION 

The X-ray emitting hot intra-cluster medium (ICM) of clus- 
ters of galaxies is known to contain a large amount of metals: 
for rich clusters between red-shift 0.3-0.4 and the present day 
the observed metallicity is about 1/3 the solar value (Mushotzky 
& Lowenstein 1997; Fukazawa et al. 1998; Allen & Fabian 
1998; Delia Ceca et al. 2000; Ettori et al. 2001), suggesting 
that a significant fraction of the ICM has been processed into 
stars already at intermediate red-shifts. 

While the origin of the metals observed in the ICM is clear 
(they are produced by supernovae), less clear is the transfer 
mechanism of these metals to the ICM. The main mechanisms 
that have been proposed for the metal enrichment in clusters 
are: enrichment of gas during the formation of the proto-cluster 
(Kauffman & Chariot 1998); ram pressure stripping of metal 
enriched gas from cluster galaxies (Gunn & Gott 1972; Toni- 
azzo & Schindler 2001); stellar winds AGN- or SN-induced 
in Early-type galaxies (Matteucci & Vettolani 1988; Renzini 
1997). 

Spatially resolved abundance measurement in galaxy clus- 
ters are of great importance because they can be used to mea- 
sure the precise amounts of metals in the ICM and to constrain 
the origin of metals both spatially and in terms of the different 
contributions of the two different type of SNe (SNII and SNIa) 
as a function of the position in the cluster. The first two satel- 
lites able to perform spatially resolved spectroscopy, ASCA and 
BeppoSAX have revealed abundance gradients in cD clusters 
(Dupke & White 2000b; De Grandi & Molendi 2001), in partic- 
ular M87/Virgo (Matsumoto et al. 1996; Guainazzi & Molendi 
2000), and variations in Si/Fe within a cluster (Finoguenov et 
al. 2000) and among clusters (Fukazawa et al. 1998). Since the 
SNIa products are iron enriched, while the SNII products are 
rich in a elements, such as O, Ne, Mg and Si, the variations in 
Si/Fe suggest that the metals in the ICM have been produced 
by a mix of the two types of SNe. The exact amount of these 
mix still remain controversial: Mushotzky et al. (1996) and 



Mushotzky & Lowenstein (1997) showed a dominance of SNII 
ejecta, while other works on ASCA data (Ishimaru & Arimoto 
1997; Fukazawa et al. 1998; Finoguenov et al. 2000; Dupke & 
White 2000a), still indicating a predominance of SNII enrich- 
ment at large radii in clusters, do not exclude that as much as 
50% of the iron in clusters come from SNIa ejecta in the inner 
part of clusters. 

M87 is the cD galaxy of the nearest cluster and its high flux 
and close location allows, using the unprecedented combina- 
tion of spectral and spatial resolution and high throughput of 
the EPIC experiment on board XMM-Newton, a detailed study 
of the ICM abundance down to the scale of the kpc. Through- 
out this paper, we assume Hq = 50kms _1 Mpc _1 , qo = 0.5 and at 
the distance of M87 1' corresponds to 5 kpc. 

Recently, when we were finishing the writing of the paper, 
we have learned of the Finoguenov et al (2001) analysis of the 
same data. Our work is complementary with theirs in the sense 
that we use 10 annular bins instead of 2, fully exploiting the 
quality of the XMM data and taking into account the multi tem- 
perature appearance of spectra in the inner regions with ade- 
quate spectral modeling. 

The outline of the paper is as follows. In section 2 we give 
information about the XMM observation and on the data prepa- 
ration. In section 3 we suggest the use of a unique set of "stan- 
dard" abundances. In section 4 we describe our spectral mod- 
eling. In section 5 we present spatially resolved measurements 
of metal abundances and abundance ratios and in section 6 we 
discuss our results. A summary of our conclusions is given in 
section 7. 

2. OBSERVATION AND DATA PREPARATION 

M87/Virgo was observed with XMM-Newton (Jansen et al. 
2001) during the PV phase with the MOS detector in Full Frame 
Mode for an effective exposure time of about 39 ks. Details on 
the observation have been published in Bohringer et al. (2001) 
and Belsole et al. (2001). We have obtained calibrated event 
files for the MOS1 and MOS2 cameras with SASv5.0. Data 
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were manually screened to remove any remaining bright pixels 
or hot column. Periods in which the background is increased 
by soft proton flares have been excluded using an intensity fil- 
ter: we rejected all events accumulated when the count rates 
exceeds 15 cts/lOOs in the [10 - 12] keV band for the two MOS 
cameras. 

We have accumulated spectra in 10 concentric annular re- 
gions centered on the emission peak extending our analysis out 
to 14 arcmin from the emission peak, thus exploiting the entire 
XMM field of view. We have removed point sources and the 
substructures which are clearly visible from the X-ray image 
(Belsole et al. 2001) except in the innermost region, where we 
have kept the nucleus and knot A, because on angular scales 
so small it is not possible to exclude completely their emis- 
sion. We prefer to fit the spectrum of this region with a model 
which includes a power law component to fit the two point like 
sources. We include only one power law component due to the 
similarity of the two sources spectra (Bohringer et al. 2001). 
The bounding radii are 0'-0.5', 0.5'- 1', l'-2', 2'-3', 3'-4', 4'-5', 
5'-7', 7'-9', 9'-ll' and ll'-14'. The analysis of the 4 central 
regions within 3 arcmin was already discussed in Molendi & 
Gastaldello (2001). 

Spectra have been accumulated for MOS 1 and MOS2 inde- 
pendently. The Lockman Hole observations have been used for 
the background. Background spectra have been accumulated 
from the same detector regions as the source spectra. 

The vignetting correction has been applied to the spectra 
rather than to the effective area, as is customary in the analysis 
of EPIC data (Arnaud et al. 2001). Spectral fits were performed 
in the 0.5-4.0 keV band. Data below 0.5 keV were excluded to 
avoid residual calibration problems in the MOS response matri- 
ces at soft energies. Data above 4 keV were excluded because 
of substantial contamination of the spectra by hotter gas emit- 
ting further out in the cluster, on the same line of sight. 

As discussed in Molendi (2001) there are cross-calibration 
uncertainties between the spectral response of the two EPIC in- 
struments, MOS and PN. In particular for what concern the soft 
energy band (0.5-1.0 keV) fitting six extra-galactic spectra for 
which no excess absorption is expected, MOS recovered the N H 
galactic value, while PN gives smaller N H by 1 - 2 x 10 2() cm" 2 . 
Thus we think that at the moment the MOS results are more 
reliable than the PN ones in this energy band, which is crucial 
for the O abundance measure. For this reason and for the better 
spectral resolutions of MOS, which is again important in deriv- 
ing the O abundance, we limit our analysis to MOS data. 

3. SOLAR ABUNDANCES 

The elemental abundances of astrophysical objects are usu- 
ally expressed by the relative values to the solar abundances. 
The so-called solar abundances can be either "meteoritic" or 
"photospheric". 

This distinction between "meteoritic" and "photospheric" so- 
lar abundances was made in the review by Anders & Grevesse 
(1989). Significant discrepancies exist between the two sets of 
abundances quoted in that paper, particularly for iron and this 
has caused in the past some controversy in the discussion of the 
results of cluster abundances (Ishimaru & Arimoto 1997; Gib- 
son et al. 1997). However recent photospheric models of the 
sun indicate that photospheric and meteoritic abundances agree 
perfectly and the community has converged toward a "standard 
solar composition" (Grevesse & Sauval 1998), with sugges- 
tions to the astrophysical community to accept this new state 



of the art (Brighenti & Mathews 1999). For the above reason, 
in this paper we shall adopt the Grevesse & Sauval (1998) val- 
ues. Since the solar abundance table used by default in XSPEC 
is based on photospheric values of Anders & Grevesse (1989), 
we have switched to a table taken from the data by Grevesse 
& Sauval (1998) by means of the XSPEC command ABUND. 
In general a simple scaling allows to switch from one set of 
abundances to the other. 

4. SPECTRAL MODELING AND PLASMA CODES 

All spectral fitting has been performed using version 1 1 .0. 1 
of the XSPEC package. 

All models discussed below include a multiplicative com- 
ponent to account for the galactic absorption on the line of 
sight of M87. The column density is always fixed at a value 
of 1.8 x 10 20 cm~ 2 , which is derived from 21cm measurements 
(Lieu et al. 1996). Leaving N H to freely vary does not improve 
the fit and does not affect the measure of the oxygen abundance, 
which could have been the more sensitive to the presence of ex- 
cess absorption. The N H value obtained is consistent within the 
errors with the 21cm value. 

The temperature profile for M87 (Bohringer et al. 2001) 
shows a small gradient for radii larger than ~ 2 arcmin and 
a rapid decrease for smaller radii. Moreover, as pointed out 
in Molendi & Pizzolato (2001) all spectra at radii larger than 
2 arcmin are characterized by being substantially isothermal 
(although the spectra of the regions between 2 and 7 arcmin 
are multi temperature spectra with a narrow temperature range 
rather than single temperature spectra), while at radii smaller 
than 2 arcmin we need models which can reproduce the broad 
temperature distribution of the inner regions. 

We therefore apply to the central regions (inside 3 arcmin) 
three different spectral models. 

A two temperature model (vmekal + vmekal in XSPEC and 
model II in Molendi & Gastaldello (2001) using the plasma 
code MEKAL (Mewe et al. 1985; Liedahl et al. 1995). This 
model has 15 free parameters: the temperature and the normal- 
ization of the two components and the abundance of O, Ne, Na, 
Mg, Al, Si, S, Ar, Ca, Fe and Ni, all expressed in solar units. 
The metal abundance of each element of the second thermal 
component is bound to be equal to the same parameter of the 
first thermal component. This model is used (e.g. Makishima 
et al. (2001) and refs. therein) as an alternative to cooling-flow 
models in fitting the central regions of galaxy clusters. 

A "fake multi-phase" model (vmekal + vmcflow in XSPEC 
and model III in Molendi & Gastaldello (2001). This model 
has 15 free parameters, as the two temperature model, because 
the maximum temperature T max is tied to the vmekal component 
temperature. As indicated in recent papers (Molendi & Pizzo- 
lato 2001; Molendi & Gastaldello 2001) this model is used to 
describe a scenario different from a multi-phase gas, for which 
it was written for: the gas is all at one temperature and the 
multi-phase appearance of the spectrum comes from projection 
of emission from many different physical radii. A more correct 
description will be given by a real deprojection of the spectrum 
(Pizzolato et al., in preparation). 

The third model is the analogue of the vmekal two tempera- 
ture model using the plasma code APEC (Smith et al. 2001). 
This model has 14 free parameters, one less than the corre- 
sponding model using vmekal because APEC misses the Na 
parameter. 

We can't adopt an APEC analogue of the fake multi-phase 
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model because the cooling flow model calculating its emission 
using APEC is still under development. Given the substantial 
agreement between 2T and fake multi-phase model (Molendi 
& Gastaldello 2001), we can regard the 2T APEC results as 
indicative also for a fake multi-phase model. 

For the spectrum accumulated in the innermost region we in- 
cluded also a power law component to model the emission of 
the nucleus and of knot A. 

For the outer regions (from 3 arcmin outwards) we apply sin- 
gle temperature models: vmekal using the MEKAL code, with 
13 free parameters and vapec using the APEC code, with 12 
free parameters. 

As pointed out by the authors of the new code, cross- 
checking is very important, since each plasma emission code 
requires choosing from a large overlapping but incomplete set 
of atomic data and the results obtained by using independent 
models allows critical comparison and evaluation of errors in 
the code and in the atomic database. 

5. RESULTS 

5.1. Abundance measurements and modeling concerns 

The X-ray emission in cluster of galaxies originates from 
the hot gas permeating the cluster potential well. The contin- 
uum emission is dominated by thermal bremsstrahlung, which 
is proportional to the square of the gas density times the cooling 
function. From the shape and the normalization of the spectrum 
we derive the gas temperature and density. In addition the X- 
ray spectra of clusters of galaxies are rich in emission lines due 
to K-shell transitions from O, Ne, Mg, Si, S, Ar and Ca and 
K- and L-shell transitions from Fe and Ni, from which we can 
measure the relative abundance of a given element. 

In Figure 1 we show the data of the 3 '-4' bin together with 
the best fitting model calculated using the MEKAL code. The 
model has been plotted nine times, each time all element abun- 
dances, except one, are set to zero. In this way the contribution 
of the various elements to the observed lines and line blends be- 
come apparent. In the energy band (0.5-4 keV) we have adopted 
for the spectral fitting, the abundance measurements based on 
K-lines for all the elements except for Fe and Ni, for which the 
measure is based on L-lines. The K-lines of O, Si, S, Ar and Ca 
are well isolated from other emission features and clearly sepa- 
rated from the continuum emission, which are the requirements 
for a robust measure of the equivalent width of the lines and 
consequently of the abundances of these elements. The Fe-L 
lines are known to be problematic, because the atomic physics 
involved is more complicated than K-shell transitions (Liedahl 
et al. 1995), but from the very good signal of XMM spectra 
and from the experience of ASCA data (Mushotzky et al. 1996; 
Hwang et al. 1997; Fukazawa et al. 1998) we can conclude that 
the Fe-L determination is reliable. Some of the stronger Fe-L 
lines due to Fe XXII and Fe XXIV are close to the K-lines of 
Ne and Mg, respectively and blending can lead to errors in the 
Ne abundance and, to a smaller extent, to the Mg abundance 
(Liedahl et al. 1995; Mushotzky et al. 1996). Also the Ni mea- 
sure is difficult due to the possible confusion of its L-lines with 
the continuum and Fe-L blend. 

5.2. Abundance profiles 

In Figure 2 we report the MOS radial abundance profiles for 
O (top panel), Si (middle panel) and Fe (bottom panel), in Fig- 
ure 3 those for Mg (top panel), Ar (middle panel) and S (bottom 
panel), in Figure 4 those for Ne (top panel), Ca (middle panel) 



and Ni (bottom panel). We note that the measurements obtained 
using the two different plasma codes agree for what concerns 
Fe, Ar and Ca; they are somewhat different for what concerns 
O, S and Ni and in complete disagreement for Mg and Ne. 

The temperature profile obtained with the two codes is 
showed in Figure 5: there are some differences in the inner 
regions, while in the outer isothermal bins there is substantially 
agreement. 

The models using the APEC code give a systematically 
worse description than the ones using MEKAL code: for the 
multi temperature models with APEC the \ 2 range from 403 
to 498 for v = 218 (216 in the central bin due to the two addi- 
tional degree of freedom of the power-law component), while 
for MEKAL models (2T and fake multi-phase give the same 
results) the \ 2 range from 323 to 382 for v = 217 (215 in 
the central bin due to the two additional degree of freedom of 
the power-law component); for single temperature models with 
APEC the \ 2 range from 546 to 886 for v = 222, while for 
MEKAL models the \ 2 range from 326 to 731 for v = 221. In 
Figure 6 we compare the residuals in the form of A% 2 between 
a 2T model using the MEKAL code and the APEC code for the 
inner bin l'-2', as well as a IT model using the MEKAL code 
and the APEC code for the "isothermal" bins 3'-4' and 1 1'-14'. 
It's evident that in the external bins the differences in the fit 
between the two codes are due to APEC over-prediction of the 
flux of Fe-L lines from high ionization states, considering the 
fact that the temperature obtained by the two codes are nearly 
coincident. For the inner regions the differences between the 
two codes are further complicated by the different temperature 
range they find for the best fit. In general, where the temper- 
ature structure is very similar, as in the innermost bin, the dif- 
ference is as in the outer bins in the high energy part of the 
Fe-L blend, while where the temperature structure is different, 
as in the l'-2' shown in Figure 6, the differences between the 
two codes are primary due to different estimates of the flux of 
He-like Si-K line. 

We therefore choose as our best abundance profiles those ob- 
tained with a 2T vmekal fit for the central regions and with a IT 
vmekal fit for the outer regions. In Table 1 we report the abun- 
dance profiles for O, Ne, Na, Mg, Al, Si, S, Ar, Ca, Fe and Ni 
obtained in this way. Abundance gradients are clearly evident 
for Fe, Si, S, Ar and Ca; O, Mg and Ni show evidence for an en- 
hancement in the central regions while only Ne is substantially 
flat. 

5.3. Comparison with Finoguenov et al. (2001 ) results 

We made a direct comparison of our results with the abun- 
dances and abundance ratios derived by Finoguenov et al 
(2001). Their values are consistent within the errors with ours 
except for the oxygen abundance, which is roughly two times 
higher in Finoguenov's analysis. To better understand the origin 
of this discrepancy we have extracted spectra from the same ra- 
dial bins, I' -3' and 8'- 14', and fitted a single temperature model 
(vmekal) in the 0.5-10 keV band, excising the 0.7-1.6 keV en- 
ergy range, in order to avoid the dependence from the Fe L- 
shell peak, as done by Finoguenov et al (2001). Our results for 
the O abundance, given in units of the solar values from An- 
ders & Grevesse (1989) to make a direct comparison with the 
results of Table 1 of Finoguenov et al (2001), are 0.32 ± 0.03 
for the l'-3' bin and 0.20 ± 0.02 for the 8'- 14' bin, to be com- 
pared with 0.535!^ and 0.386^;^. Also interesting are the 
results for the Ni abundance: 2.55 ±0.33 for the l'-3' bin and 
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2.34 ±0.26 for the 8'- 14' bin, to be compared with 2.573^;^ 
and 0.800^ 800- We also note that a 2T modeling of the inner 
bin I' -3' gives a statistically better fit over a single temperature 
model ( x 2 /d.o.f of 919/440 respect to 1083/442, using the 0.5- 
10 keV band) breaking down the assumption of near isother- 
mality. 

We also cross compare the K- and L-shell results for Ni and 
Fe in our analysis, performing the spectral fits for all the ra- 
dial bins in the 0.5-10 keV band, but excising the 0.7-1.6 keV 
energy range. The derived abundances for the two elements 
are consistent within ler, although the K-shell Ni abundance is 
20% higher than the L-shell abundance in the bins fitted with a 
single temperature model. It should be borne in mind that par- 
ticularly for the outer bins the K-shell Ni abundance measure is 
very sensitive to the background estimate. 

5.4. Abundance ratios and SNIa Fe mass fraction 

From the abundance measurements we obtain the abundance 
ratios between all the elements relative to Fe, normalized to the 
solar value. They are shown in Figure 7, Figure 8, Figure 9 and 
in Table 2, together with the abundance ratios obtained by mod- 
els of supernovae taken by Nomoto et al. (1997) and rescaled to 
the solar abundances reported in Grevesse & Sauval (1998). We 
use those abundance ratios to estimate the relative contributions 
of SNIa and SNII to the metal enrichment of the intra-cluster 
gas. Such estimates are complicated by uncertainties both in 
the observations and in the theoretical yields. Our approach is 
to use the complete set of ratios trying to find the best fit of the 
function 

xXq/Fcq/ observed V Xq/Fcq J SNIa \ Xq/Fcq ) SMI 

(1) 

where ( v X {t e J is the measured abundance ratio of 

\ x e/ Fe @ J observed 

the X element to Fe, given in solar units, ( Y X iC e ] and 

\ X o/ Fe eJ SNIa 

( v X i^ e I are the theoretical abundance ratio by the two 

types of supernovae, also given in solar units and / is directly 
the SNIa Fe mass fraction. The result of the simultaneous fit 
of the eight ratios is presented as circles in Figure 10, using 
the SNII model by Nomoto et al. (1997) and W7, WDD1 and 
WDD2 models for SNIa respectively for the three panels. Due 
to the large uncertainties in the yields for the SNII model, the 
results are strongly SNII model dependent. For comparison we 
use the range of SNII yields calculated by Gibson et al. (1997), 
also listed in Table 2, which involve only ratios for O, Ne, Mg, 
Si and S and as before finding the best fit for eq.(l). The re- 
sults are shown as squares for the lower end of the range and 
triangles for the upper part of the range. 

The best fits are obtained in the inner bins with a combina- 
tion of the WDD2 model for SNIa and the Nomoto model for 
SNII, with reduced x 2 which ranges from 2 in the inner bin up 
to 10 in the 4'-5' bin. In the outer bins the fit is slightly better 
(reduced x 2 of 8-10 instead of 12-13) with a combination of 
the W7 model for SNIa and the Nomoto model for SNII. This 
is shown in Figure 1 1 where the fits with the three SNIa mod- 
els together with the Nomoto SNII model are reported for the 
0.5'- 1' bin and for the 11 '-14' bin. It's clear from the inspection 
of the residuals in terms of A% 2 that the combination of W7 
model and Nomoto SNII model fails in the inner bins because 
it predicts a higher Ni/Fe and O/Fe ratio, compared to the de- 



layed detonation models. The situation is the opposite for the 
outer bins where a higher Ni/Fe and a lower S/Fe favors the W7 
model. However the preference for the W7 model in combina- 
tion with the Nomoto SNII model in the outer bins is strongly 
dependent on the Ni/Fe ratio: if we exclude it from the fit the 
WDD2 model provides the better fit to the data. 

6. DISCUSSION 

The model emerging from the ASCA and BeppoSAX data for 
the explanation of abundance gradients in galaxy clusters was 
that of a homogeneous enrichment by SNII, the main source 
of a elements, maybe in the form of strong galactic winds in 
the proto-cluster phase and the central increase in the heavy el- 
ement distribution due to an enhanced contribution by SNIa, 
strongly related to the presence of a cD galaxy (Fukazawa et 
al. 1998; Dupke & White 2000b; Finoguenov et al. 2000; De 
Grandi & Molendi 2001; Makishima et al. 2001). As a text- 
book example we can consider the case of A496 observed by 
XMM-Newton (Tamura et al. 2001). The O-Ne-Mg abundance 
is radially constant over the cluster, while the excess of heavy 
elements as Fe, Ar, Ca and Ni in the core is consistent with the 
assumption that the metal excess is solely produced by SNIa 
in the cD galaxy. The crucial ratio for the discrimination of 
the enrichment by the two types of supernovae, O/Fe, is then 
decreasing towards the center. 

The XMM results for M87/Virgo question this picture. They 
confirm and improve the accuracy of the measure of heavy ele- 
ments gradients previously found by ASCA and BeppoSAX, but 
they also show a a statistically significant enhancement of a 
elements O and Mg in the core. If we consider the inner 9 ar- 
cmin the ratio O/Fe is constant with x 2 = 6-9 for 7 d.of. and 
adding a linear component does not improve the fit (x 2 = 5.4 
for 6 d.o.f). These facts points toward an increase in contribu- 
tion also of SNII, since O is basically produced only by this 
kind of supernovae. Although there is little or no evidence of 
current star formation in the core of M87, the O excess could 
be related to a recent past episode of star formation triggered 
by the passage of the radio jet, as we see in cD galaxies with 
a radio source (A1795 cD: van Breugel et al. 1984 and A2597 
cD: Koekemoer et al. 1999), nearby (Cen A: Graham 1998) and 
distant radio galaxies (van Breugel et al. 1985; van Breugel & 
Dey 1993; Bicknell et al. 2000) (for a comprehensive discus- 
sion see McNamara 1999). 

To put the above idea quantitatively, waiting for a true depro- 
jection of our data, we use previous ROSAT estimate of the 
deprojected electron density in the center of the Virgo clus- 
ter (Nulsen & Bohringer 1995) to calculate the excess mass 
of oxygen. To estimate the excess abundance we fit the in- 
ner bins, where we see the stronger increase in the O abun- 
dance, with a constant obtaining an abundance of 0.32, while 
for the outer bins we obtain an abundance of 0.21 so the ex- 
cess is 0.11. Then we estimate the oxygen mass, Mo to be 
M = A o yo !0 Z e xcess M H , where M H = 0.82» e f {R out -R in f, R in 
and R ou , being the bounding radii in kpc of the bins, A = 16 and 
yo,Q = 6.76 x 10" 5 from Grevesse & Sauval (1998). We obtain 
a rough estimate of 10 7 M Q which, assuming that about 80 M Q 
of star formation are required to generate a SNII (Thomas & 
Fabian 1990) and Nomoto SNII oxygen yield, requires a cumu- 
lative star formation of 4-5 x 10 7 M Q . This star formation is 
in agreement with a burst mode of star formation (< 10 7 yr) at 
rates of ~ 10-40M Q yr _1 , as it is observed in the CD galaxies 
of A1795 and A2597 (McNamara 1999). 
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For what concern heavy element gradients and the contribu- 
tion of SNIa, M87 data suggests an agreement with delayed 
detonation models (in particular for the inner bins), as stressed 
by Finoguenov et al (2001), in contrast with the preference of 
W7 model set by the high Ni/Fe ratios found by Dupke & White 
(2000a). This is particular evident if we consider the S/Fe ra- 
tio in Figure 7. If we consider the set of theoretical values for 
SNII and W7 SNIa models the behavior of these ratio would 
indicate an increasing contribution by SNIa going outward to 
the center. We recover the correct behavior if we choose the 
WDD1 yield and we reduce the S SNII yield of Nomoto et al. 
(1997) by a factor of two to three, as was already indicated by 
ASCA data (Dupke & White 2000a). We caution however that 
this is a substantial contribution larger than that allowed by the 
SNII models choosed by Gibson et al. (1997). The use of de- 
layed detonation model for SNIa could also explain the over 
abundance of S and also Si (respect to the W7 model) found by 
Tamura et al. (2001) for the core of A496. 
With increasing radius the W7 model gives a better fit to the 
data respect to delayed detonation models. This fact indicates a 
SNIa abundance pattern change with radius and could be taken 
as an independent X-ray confirmation of the conclusions of 
Hatano et al. (2000) on the optical spectroscopic diversity of 
SNIa, as suggested by Finoguenov et al (2001). However we 
stress that the preference of W7 over WDD models in the outer 
bins is entirely due to the Ni/Fe ratio which could be affected 
by systematic uncertainties, as discussed in section 5.1 and 5.3. 

In Figure 10, we show the relative importance of SNIa, mea- 
sured by the Fe mass fraction provided by this kind of super- 
novae. This is substantially constant through the 14 arcmin an- 
alyzed. This fraction is considerable and ranges between the 
50% and 80% and it depends only slightly on the SNIa model 
used. Instead the uncertainties involved in using different SNII 
models are large and a definitive answer cannot be reached until 
further convergence of SNII models is achieved. 



7. SUMMARY 

We have performed a spatially resolved measurement of the 
element abundances in M87, the Virgo cluster cD galaxy. The 
main conclusion of our work are: 

• the APEC code gives a systematically worse description 
than the MEKAL code in modeling M87 spectra; 

• we confirm the increase of Fe and other heavy elements 
towards the core indicating that the SNIa contribution 
increases; 

• the increase in O abundance and a constant O/Fe ratio 
in the inner 9 arcmin indicates an increase also in SNII 
ejecta possibly from star-burst in the recent past; 

• Si/Fe and S/Fe profiles favor WDD models over W7, 
also requiring substantial reduction of the SNII yield of 
S; 

• the indication of a change of the SNIa abundance pat- 
tern, provided by a preference of the W7 model over 
delayed detonation models in the outer bins, is entirely 
due to the Ni/Fe ratio. Since the Ni measurement is dif- 
ficult and uncertain, this indication should be taken with 
some caution. 
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FIG. 1 . — Data of the 3'-4' bin and lines of the various elements, obtained by setting all element abundances to zero except the one of interest, calculated using 
the MEKAL code. Top Panel: the lines of O, Fe, Ne, Ni, Mg and Si used for the analysis in the 0.5-2 keV energy range. Bottom Panel: the lines of Si, S, Ar and 
Ca used for the analysis in the 2-4 keV energy range. 
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FIG. 2. — MOS abundance profiles for O, Si and Fe. Uncertainties are at the 68% level for one interesting parameter(Ax 2 = 1)- Diamonds represent measurements 
with 2T model using APEC code, asterisks those with 2T model using MEKAL code and squares those with fake multi-phase model using MEKAL code. Full 
circles indicate measurements with IT model using MEKAL code while empty triangles those with IT model using APEC code. 
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Table 1 

Abundance profiles obtained by a 2T vmekal modeling for the central regions and a IT vmekal modeling for the 

OUTER REGIONS. ALL ABUNDANCES REFERRED TO SOLAR UNITS AS GIVEN IN GREVESSE & SAUVAL (1998). ALL ERRORS QUOTED ARE 

AT THE 68% LEVEL FOR ONE INTERESTING PARAMETER (Ay 2 = 1) 
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Table 2 

Abundance ratios from our results and from SNe models. The abundances relative to Fe normalized to the solar 
value, (x/fe)/(x /fe ), where x and fe are number density of the element and fe, are shown. the errors associated 

with the observed abundance ratios are the propagated la errors. 



0'-0.5' 
0.5'-l' 
l'-2' 
2'-3' 
3'-4' 
4'-5' 
5'-7' 

r-9' 

9'-W 
11 '-14' 



O 

30 +0 06 

u - JU -0.04 

39+ 03 

T,A +0M 
u - J ^-0.03 

0.39+2-2J 



0.31 
0.40 



+o:o4 
: 8:o4 

-0.04 



n qo+o!o4 

041 +0.04 
u>fi -0.04 

45+ 05 

u -^°-0.06 

52 +0 - 05 
0.031 



Ne 
29 +ai3 
35 +01 ° 

n 4A+O.IO 
U - W -0.07 

45+ 07 

u -^°-0.07 

69 +008 

74+0 07 
u - '^-0.06 
57+0.06 
u ""-0.06 

61 +0 08 

u - J -'-0.08 



Mg 
62+ 009 
56+ 004 
60+ 004 
62+ 004 

u - oz -0.04 

64+ 004 

u -°^-0.06 

0.743* 
0.88^ 

0.89^852 

56+ 008 

U - JO -0.07 

52+ 009 



Si 

OQ+0.12 

uu -o.io 
6 +o.04 

uo -0.05 
i -7+0.04 
iz -0.05 
04+0.04 
UH -0.04 
i o+0.04 

lol 

-0.04 
fi +0.03 
UD -0.03 
00+0-05 
uu -0.04 

ns+o.04 

U8 -0.04 

2 6 +o.05 
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n oa+0.10 
U -<"-0.08 

76+ 03 
u - '"-0.04 
2+0.03 
-0.03 

75+ 03 

u - ' J -0.03 

74+ 03 

J -0.03 

63+ 03 
u - D;) -0.03 

o sn+o- 03 

4S+ 04 

u -^ J -0.04 
j +0.04 
-0.05 



0.73! 



0.65 + 



0.33! 



Ar 

36+0.21 
■ JO -0.20 

51 +o.io 



48- 1 



-0.15 
+0.12 
-0.12 
cc+0.12 
JJ -0.12 
cn+0.15 
J ^-0.14 
6O+0.15 
oy -0.16 
cc+0.13 
JJ -0.13 
07+0.15 
u/ -0.16 
Ti+0.19 
Ji -0.21 
c 4 +0.22 
"^+-0.24 



Ca 

1 19+O.24 
1 - lz -0.23 
1 90+O.I6 

60+014 

• DU -0.14 
_ 37 +0.13 



1.82 + 



1. 

1 69+" : ^ 

-0.16 

1 97+0.16 
1 - z/ -0.17 
1 C4+O.2O 
i - J ^-0.21 

1 50 +0 ' 23 

1 ■-' u -0.25 



Ni 

1 41+0.25 

1 61 +0 ' 15 
i - oi -0.12 

D+0.11 

s -0.10 

1 60+ 012 

0.12 

.8:14 
-0.12 
1 74+0.12 



1.48! 



1.51! 
1.69! 



-0.11 

9 17+O.n 
z - 1 '-0.14 
9 41+0.21 
z -^ 1 -0.16 
9 9O+0.22 
Z - /a -0.18 

3.23 



W7 a 



0.004 



0.022 



0.36 



0.30 



0.42 



0.33 



WDDl fl 



0.030 



0.002 



0.036 



1.15 



1.01 



1.56 



1.48 



1.12 



WDD2 fl 



0.016 



0.001 



0.013 



0.69 



0.62 



0.98 



0.97 



0.95 



SNIP 



3.23 



1.90 



2.58 



2.37 



1.17 



1.72 



1.09 



1.27 



SNIF 



1.25-2.80 



0.96-1.93 



0.96-1.89 



1.81-2.54 



0.85-2.08 



a Different models of SNIa taken by Nomoto et al. (1997). 
* Yields of SNII taken by Nomoto et al. (1997). 

c Gibson et al. (1997) who choose a representative sample of SNII yields in literature. 



